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Re CODE NUMBERS 


A system of coding data points to reflect in a single 
multi-digit number, tube geometry, angle of inclination and 
sequential number of the test run was developed. The run code 
momoer is a four digit number; the first number indicating 
moe geometry, the second indicating tne angle of inclination, 
ame. the third and fourth numbers indicating the sequential 


moumpober of the run. 


Tube geometry codes: 


pode . Desert ion 
ak Staqgered tubes, Sieoncl =o Dy (Smith-Kiss and 
Bond-wallin data) 
2 Ln=1ine Cubes 50 t-e oa Dy 
3 Staggered tubes vs? = 125 Dy cea 2) D, 
4 In-line tubes “on =) iy s Dy SL = 2.0 D 
Angle codes: 
Code Desceoelon 
5 30 degree angle of inclination 
4 45 degree angle of inclination 
6 60 degree angle of inclination 
y 75> Geqree. anGle) Of Vinecli netic 
2, 90 degree angle of inclination 


Example: 1704 


EUDeG Geometry mil, of wawo ls 9), staguered tubes. 


i? 


seventy five degree angle of inclination, 
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While adequate empirically derived design information is 
available for heat transfer equipment where the tube banks ere 
Meeallel to or normal to the flow of the fluid being heated exe 
cooled, little information exists for the case where the tubes 
are inclined at some angle to the fluid flow. The cése of 
mabes at an angle to fluid flow is not uncommon perticulerly 
Mmeemarine boller convection banks. It is the purpose of this 
Mhesis to investigate the effect of variation in flow angle on 
Menvection heat transfer and frictional resistance in a tube 
Meteor stig 2 moael Scale apparatus, and to Gevelop ec iatencs 


Tomer Will be useful to the designer of Heat transfer equipmenmas 


The relationship originully developed by Grimiscn cal and 
Meer modified by the transport data obtained by Justi and 
Leuder[ 2] will be us'd as the basic convection heat trénsfer 
Meraicionshnip where the tube banks are inclined at minety Gages 
Meee fluid flow. A feccor Bo will be developed which, wher 
incluued in the moaificd version [2 Jof this ¢«.u.tic. will give 


ie che fort: 





one See oh Ge0U BN ee 
=< Umcogc eM Peta ye Ps . / 
K 2h od Oo ie Ch A os i Ms 
ie mac Lei ae Wild be Ob Gol ned {renaGieese CULVes. send wien Gem 
used will give a modified value of Nusselt Number np, tor Ube 





tubes inclined at an angle ©@ to fluid flow. ‘ 


CONCUL sent ly, Swi ome Slew tnVe sel G ee met he wet Lec ee: 
Variatlon/in flow anole cpr convec tion, Wat trensier cie eLiece mam 
frictional resistance as expressed Byrrhe Fanning type frictrem 


factor will be determined. 








PNumbers in brackecs Geciagqnete referene@tS listed atache ena ie. texas 


—_— 
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The availability of electronic computers makes possible 
more accurate reduction and analysis cof experimental data 
than has been possible before. In addition, they permit a 
@reater number of calculations to be carried out in a short 
time. Thus a systematic series of geometries wiil be tested 


by varying the longitudinal tube pitch and tube arrangement. 
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THEO 


Effect on Convection Heat Transfer 


The transfer of heat by a fluid passing through a tube 
bank is by the mechanism of forced convection in turbulent flow. 
The quantity of heat transferred through the stagnant fluid film 
Or thermal boundary layer contiguous to the tube is a function 


of the following quantities: 


a) a characteristic length, taken in this development as 
the outside diameter of the tubes, Dis 

b) the coefficient of convection heat transfer, Des between 
the fluid and the tube surface. 


@) “the “thermal ‘conductivity i. .con Pthe cr Eid, 
d) the specific heat, Sa of the fluid at constant pressure. 


ej thes aynemice Viscosity .<- wor sthe f£iuad, 


l- 


£f) the mass flow of the fluid, G. 


The convection heat transfer coefficient may be expressed 
\ 


as a function of the other quantities: 


~ hs» £ (Dt, K, cp,#, G) 


The choice of a characteristic length and of a Cemperarcume 
at which to evaluate those temperatures dependent physical 
properties is open to some discussion. when fluid flows across 
a single tube the outside diameter of the tube is the most 


logical choice, Nowever in ja tupe=bank the flow 1s veriectea py 





ee 
the clearance between tubes and other definitions of character- 


metic length have been suggested. 


LJ 


To preserve consistency with Pierson, ememte! 4! une aueua 
and because other definitions provide no proven edvantage the 
Outsiage diameter of the tube will be used. Similarly various 
Mefinitions of the temperature at which to evdluate physical 


properties have been suggested. 


That which is consistent with the work menticned above 
‘emvicl 25 Most 1Ogreal For Loreced Convection ical teanstes 
is the temperature of the fluid film. As there is a variation 
moe tua film temperature through its thickness the film 
temperature is defined as the arithmetic mean between the bulk 
Gemperature and the tube surface temperature. Likewise the 
bulk temperature varies through the tube bank and is defined 


as the arithmetic mean of the inlet ana exhaust temperatures. 


The exact mechanism of forced ccnvection heat transfer 1s 
moe fully understoud so that while analytical methods can 
sometimes be used, semi-empirical or empirical relaticnships 
are more reliable. The uevelopment of empirical relationships 
is asSisted by the derivation of dimensionless groups of the 


physical qualities discussed above. 


These groups may be derived by the methods of dimensional 
analysis. The applicable dimensions are: 
ite fOr Seng an 
T for time 
Mion OF Mae S 
metOr Heat energy 


@ for temperature 
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AS heat 1s equivalent to mechanical work it can be 


c noe The 


replaced by ~force) x (resistance) or H « ML 
dimensions of the remainder of the physicéel properties may be 


expressed as: 


D, = L 
h, = mr >p7t 
k = mLT7?¢7t 
es : om ape 
fu = mitt 


ef= mE ep7 


Meee Buckingham Pi Theorem 1s applied, the following matrix 


is developed: 


Dy Ne K os fi G 
M O 1 1 O i 1 
iE 1 0 i 2 = =. 
T 0 ~3 =3 =O =) =1 
“y C = =i -1 O 0 


Prom Chis matrix three dimensionless relations whicn mave 


proved uséful in heat transfer work are uerived. These are: 


eae = MT7"O7* 1 = Nusselt number, 


k 
ut~-¢7+ p 


Ween may physically De Toiecerpreted aS rine ate jo, enc 
Gmaracteristic Lengen- Oa Stationary. rete ae nove tine 


heat at the same rate under the same temperature uliference; 











tt — 


19 
Sea musa = Reynolds number, 
mu7 trot 


which may be interpreted as the ratio of inertia force to 


wrsccous force and: 


1 


1 = Prandtl Number 


“pM ann 2d-l 
mir >g7t 


which is the ratio of two constants of molecular transportation. 


ML” 





The relationship is usually written in terms of power functions 


h_-D, an opt? cM b 
k fA Se 


where the coefficient A and the exponents a and b are constants. 


as: 





By correlation of empirical data where the fluica flow is 


normal to the tube bank the above equation has been shown to be: 





h -D 
oie " 0.60 Ie 
7 =a lO) Paras) 2 FF, GD. cp M 
7 k 
where Py and F, are empirical correlations known as arrangement 


Maetor and G@épth factor respectively. Ine are rangeneni ae. decem 
is a function of tube arranyerent, staggered or in-line, 
longitudinal and transverse spacing ana Reynolas Number. The 
aepth factor 15 e@ finction of the “number “of rows aneene soup 
bank. The deoth factor determined by Grimison becomes unity 


FOr LEN or more rews, 
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C®a burr”! nas Showa rthat, tor fluid flow parallel. eorvthe 


tube surface; 





Peanem=) 01025 oe Sel 0 
; Car 
K 
Where dD. is the hydraulic diameter. In order to be able to 


make a direct comparison with flow normal to the tubes the 
above expression may be written in terms of tube diameter. 


h,-D 


c 0.80 
x SOUS (=) 


0.20 1/2 
(ot) C8") 
Dy k 


For cases where the fluid flow is neither normal nor 





in 


parallel to the tube surface, but inclined at some angle theta 
€O it, the rate of heat transfer as indicated by the value of 
Nusselt Number should lie somewhere between that for normal 
flow or cross flow, and parallel flow. It should be possible 
Mmemenxplress the effect of the inclined flow on forced convecticgn 


heat transfer in one of two ways either as: 


1) a modification to the exponent of Reynolas Number, i.e.;3 


h__D 0.60 £(89) 
ies = 0,.2924F-. GD Cc 
== A “ S| ee) 


where £ (0) represents a function of the angle of inclination, or 


elm a5 an acda tional mulciplying.tacter. 


y h _D OPS 6 Ly ee 
vm k 


Se 





Ak 
An analytic prediction of the effect of flow angle on 
convection heat transfer would be at best uncertain. As the 
Beriier work with which this work is to be consistent is an 
empirical correlation, no theoretical prediction will be 
attempted. Instead, data derived from the test apparatus 
will be correlated witn the empirical equation to produce a 


feete ying factor which is a function of the angle of inclination 


ELTECeE PON UL riCr tonal shesi stance 


The frictional resistance of a fluid flowing across a 
tube bank is customarily treated by an equation .f the general 
Fanning type wherein the pressure drop is treated as being 
caused by flow through a serics of orifices formed by successive 
Matfer restrictions in the flow path. By methoas of dimensional 
analysis similar to those described above for convection heat 
transfer a dimensionless friction factor, f£f, may be derived 


which is a function of Reynolus Number only: 


f= 3) 





AA 


where the coefficient C ana tne exponent X are constants. 
This equation may be modified to permit the use of convenient 


Guantities without reducing its soundness: 


£f = 10,84 e208 


NG* 





OLE 


where: 
eo a COnStanht whicneincludes the effect of the 


deecel Cr atlLOnmge. “Gravity In ft... per hrs per hrs 
end weead "GOonversi On Lactor from pounds per 


Sseulane f£OOt EOVINGNeS Of Water 


=salr GCensi cy 


DP eS Oressuberemouracross Ene Cube voankeinotmenes 108 
water 
N SwMUNDerN Oe inciOEetadechtitecal resturteti Ona, 


aid the case of a tube bank N is equivalent to 


| the NuUMDer SLT rows wor tubes : 


G = mass fdow 


The question of what temperature should be used to determine 
Miemoceperties of the fluia again arises. The temperature has 
emeeroniticant effect on gas density and viscosity which directly 
SaeeoGit tie friction factor and the corresponding Reynolds Numbems 
Utilization of the mean bulk temperature in the bank has been 
Shown to be unsatisfactory by previous investigators. Empirical 
aevelopments led to the establishment of the convention that 
the physical properties be evaluated at a temperature equal 
femene Cube Surface temperature minus eight tenths the log mea 
Bemberdture ditrerence for staggered arrangements or equdl te 
tube surface temperature minus nine tenths the log mean 
temperature difference for in-line arrangements. The minus 
Sign applies where heat is transferred from hated tubes to 


eeeGooler fluid, 





es 


ewe tigiemvi Li) UNG Convection heat transter tite 
Seedtction Of variation of friction factor with flow angle 
Maes not tend Ltself to analytic prediction. Accordingly no 


agetempt at analytic prediction will be made. 


The choice of a purely empirical approach while not as 
intellectually stimulating as a theoretical development has 
the distinct advantage of producing results which are 
reproduceable and can be used with some confidence in the design 


of heat transfer equipment. 
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EXPERIMENTAL TECHNIQUE 


The technique followed in the investigation is that of 
forcing air through a bank of heated tubes which are arranged 
in various geometries and inclined at various angles to the 
air flow in order to measure changes in parameters indicative 
Srmiciie quantities of heat transferred and the frictional 
resistance in the bank. Data is obtained for three geometries, 
all with a pitch normal to the air flow of 1.50 tube diameters. 
Pitches in the direction of flow are 1.50 and 2.00 tube diameters 
and the tubes are arranged in staggered and in-line configura- 
tions. For each geometry, measurements are taken with the 
tube bank at ang! esvor. JO; 2o0- Aopand 30 degrees to the air 
flow. The angle is varied by inclining the tube bank ina 


test section of constant cross section and length. 


The tubes are internally heated by steam, the inlet pressure 
being controlled so that the tube surface temperature is very 
Meeeeele degrees Farenheit. Ambient. air 1s heated Dy being 
forced through the tube bank by a steam turbine driven blower. 
Variations in the mass flow of air are obtained by varying the 
speed of the supply blower. The mass flow of air is measured 
by means of a square edged orifice and appropriate manometers. 
Temperature sensors in the air duct upstream and downstream 
from the tube bank are used to determine inlet and outlet 
temperatures. Tube surface temperatures are determined by 
sensors located on the tube surfaces tnroughout the bank, 


penerally sony chev upstreamn sector ver the tube. ube: surtace 





eo 


temperatures used in the computations are then the arithmetic 
mean of the measured tube temperatures. The principle parameter 
for determining the effect of frictional resistance is the 
pressure drop across the tube bank measured by manometers. 
Transport data is evaluated at appropriate temperatures with a 
humidity correction obtained by standard psychometric techniques 


m@eeired to the specific heat. 


In addition to the data obtained from the three geometries 
tested, the data determined by previous investigetors has also 
been included. All this data has been placed through the 
computation process and analyzed for correlation of both the 


previous work and the extended investigation. 








PROG@eo Wem | OK REDUCING EXPERIMENTAL DATA 
Ticmemrw@eseror the reduction of experimental data is the 
calculation of Nusselt Numbers, Reynolds Numbers and friction 
Mero, This process may be divided into four parts: 
1) calculation of weight rate of flow through the apparatus 
2) cCaleuiation of Nusselt Number and the corresponding 
Reynolds Number 
Bo eecalcilacioOn of friction factor and the corresponding 
Reynolds Number 
A) calculation of Nusselt Numbers by the mcdified Grimiso@ 
equation ana the Colburn equation. 


This section develors in detail the calculétion proceuure. 


1) Weight rate of flow: 
The weight rate of flow calculation is perhaps the mos@ 


Setolex of those performcd. The equation is: 


WAO = A, CB) eg (P, - Po) 


Wwnere: 
WAOQ = weight rate of flow 
AK = flow coefficient 
YoI=Pexoansilon Tractor 
Gos Graviteticna: constant 
f= Sin ener Ly 
[Py-P>} = pressure differential across the orifice 


The metnod used £6 obtain weigqnte rate of flow from orituee 
dimensions and pressure differential is a standard engineering 
technique and hus been developed in detail by both Bond anu 


ae (7 
conn ere Kiss and Smith; ee Wide NO ele te Deat oa ier. 


_ 





2) 
BY GoOroing Constants and applying proper Unit conversics 


Pactors the equation reduces to: 


6) (0) 224 ES (IDO)* (AK) Y J DHF 


WAOQ = 
Pier e: 
MOOS ="Critace meter Gdiameter win inches 
DHF = pressure drop across the orifice meter in inches of “cam 


The expansion factor Y is a function of diameter ratio, 
meet Of differential pressure to absolute static pressure and 
the ratio of specific heats. As shown in Appenaix D the factor 
¥Y can be calculated for each data point by applying an equation 
of the form 


weet = C (DHF) 
Py ana 13.596PB 


The flow coefficient, AK, is a function of orifice diameter 
M@eemeter ratio, and pipe Reynolds Number. The aqerivation of 
Meemexoressions for flow coefficient in Appendix D gives 
Pevacions of the form: 

ANG C5 + C. 

Y ANRE- 
However Reynolas Numover in the pipe is given by the 


Sore sslon: 


ANRE = VW (IDO) (AK) C /DHF 


where: 


ANRE 


Pipe Reynolds Number 


V BxpanS LOM = cGgiEe j 


TDO Ortfice:- Diameter 


H 





28 
oo) myname boCOSity fOr air, the expression for which 


was developed in Appendix D and is #4 


Dee (cionit 2 12.6) elon 
AK = flow coefficient 
eee wer eonstane 
DHF = differential pressure across the thin plate orifice 


meter. 


It can be seen that flow coefficient and Reynolds Number 
are mutually dependent and an iterative process must be used 
to determine satisfactory values for these two quantities. 
Mats Can be done by first calculating intermediate quantities 
composed of those factors other than AK which appear in the 
Reynolds Number expression specifically: 


Vi= Beco? y 107 
O.16TI + 1ll2.t 


and: 


Yue Jl = Ci (DHF) 


Py + 13.596PB 


A trial value of AK depending upon orifice diameter may be 
assumed and a corresponding value of Reynolds Number may be 
@aleulated : 

ANRE = V x Y¥ x AK 

This value of Reynolds Number may be used to calculate a 
second value of flow coefficient which in turn may be used to 
compute a new Reynolds Number. This process can continue 
until the comparison of the new value of AK with tne previously 


pag 








ao 


calculcted value shows that the increase in accuracy docs not 
PeGeeant repetition. This final value of flow coefficient 
Sep then be used to celculate weight rate of flow of air 


emcough the apperetus. 
Pew ncseno76) (ibOyey (AK) <(x) 


When the weight rate of flow has been determined certain 
Other values may be daetermined. These include: the average 


temperature in the tube bank in degrees Rankine: 


Ca 


and the specific heat of wet air by the method aetailed in 


ate 459.09 


Appendix D, the log mean temperature difference: 
hisilD = te ee oS 


LUG, Lif L 
A ae Es 


ana film tenpercture: 


ae ee \ a 
one 
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when these quantities have been computed everytning 
mecessSary for the second major calculation Step, the computatzem 
cf Nusselt Number has been completed. The computation of 


Nusselt Number may then be undertaken. 


2) Nusselt Number: 
Nusselt Number may be defined as 


Nusselt Number = Q 


———— . bt 
S (XLMTD) K 
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where: 
OF= total Meat transferred 
Dt = tube diameter 
S = total heating surface 
mion4tD = log méan temperature difference 


Pes eENehMmal CONGUGLIVICY Of aly 


The relation for Nusselt Number may then be expressed as: 


ANNU = WAO (CPW) (TE = TI) 


Cl ure) (CG, 000074 40, 01344) 
where: 
eae waO (CPW) CTE=-TL) 


m= 0.00002Z4TF + 0.01344 


Biyoo= effective length L x IT 


Pyeueliiizing the interior dimension of the box, minimum possiodle 
moe ebength, and combining constants the equation reduces to: 

ANNU = wAO (CPW) (TE -— TI) 

Peob ihe OO Ieee Die Gee. 

For flow in which the tubes are inclined the surface area is 
increased. The surface area at any angle of inclinction is 
equal to the surface area at ninety degrees divided py the 
Sine of the angle of inclination, thus for angles other than 
Pmeuy degrees the expression for Nussclt Number beccmes: 


BNNU = WAO (CPW) (TE - TI) SINE 8 
XEMID COsOrPset. == 6.4679) 


= ANNU (SINE ©) 





ot 


Mees 15 the proper value of Nusselt Number for inclined flow 
and is the one result of the second part of the calculation. 
meets Now necessary to calculate the Reynolds Number through 


the test section. Reynolds Number may be expressed as: 


BNRE = WAO (Dt) , 


A sw 
where: 
A = free flow area through the test section, a projected 
area which remains constant though the angle is varied 
moe Oynemic ViSGOSlty Of air = (.16TF x 112.6) x oe 


The equation reduces to: 


BNRE = WAO x 10° 


CO. O46 25601 32.252) te 
with these two results, Nusselt Number and Reynolds Number the 


Second part of the computation 1S complete. 


See briction factor 
For the frictional resistance calculation some quantities 


which have not been previcusly calculated are required. 


KR new definition of film temperature is necessary. This 
definition is a function of tube temperature and log mean 
temperature difference and may be expressed in the form: 


Ver PT = CeeCxUMrD) 


5 
where Ce is either 0.800 or 0.900 depending upon whether the 


tube geometry is staggered or in-line. A Reynolds Number for 





a2 
woe viel) Clemente t ron Tractor may be calculated Using the 
vit Ue of film temperature just computed: 
CNRE = WAQ x 10° 
| COR O46256 TFF + 32.552). 
The friction factor is expressed as: 


FRIFAC = 10.84 (RHO) pp 10° 


m (WADE 
A 
The density of air is calculated using the equation: 


aos = RT in the form 
RHO 


fo = (Py + 13-596PB) 5.1815 Yi 
93-349 TRF = 459,69 
When RHO has been obtained the friction factor may be calculated: 


FRIFAC = 1,084 (RHO) (DP) 10° 


WAo® 
penlcys- 
4) Nusselt Number by the Grimison Equation 
When the friction factor has been obtained, the 
computation continues to calculate the Nusseit Number by the 


modified Grimison equation. 


The Grimison equation is expressed as: 


TNNU = 0.292 (ARGFAC) F 


0.60 iis 
d (BNRE) (cenit 


k 


Po-ENe CeStrcect lone su cen rows, dGepuene depUchn factomu. 
P42 2S unity and can be eliminated from the expression. The 


Mantle res Reynolds Number, BNRE, and Specifie heat, CPW are 





Sy. 


available from previous work, and analytic expressions in 

terms of temperature exist for dynamic viscosity and thermal 
conductivity. It has been shcwn in Appendix D that the 
arrangement factor may be expressed as a function of Reynolds 
Number for specific geometries. With proper suystitutions the 
expression for Nusselt Number beccmes: 

ers) 


TNNU = 0.292 ARKGFAC (BNRE) °° er (UG. lor +217.6) Gime) 


sCU002Z4ir + 0201344 


Jhen this value of Nusselt Number has been determined, a 
factor FTHLTA is defined by forming the quotient of the 
experimental Nusselt Numver BNNU, and the Grimison Nusselt 
Number, TNNU: 

FTHETA = BNNU 

TNNU. 

wnen FTHETA has been computed, the computation of Nusselt 
Number by the Colburn expression is undertaken. The similarity 
an the Colburn expression and that for tne Grimison equation 
MefesecO Simplifying the calculation by calculating the Colbuge 
Nusselt Number using the previously calculated Grimison 
Nusselt Number: 


0.80 


CNNU = .023 (BNKE) (DIAPAC) (4 


ia 


cs 
mere DIAFAC 1S a conversion factor tne ture diameler to 
hydraulic diameter. 


PVN sae SO (BNRE) 2? °8( ARGFAC) (spent t 
K 


ee 
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bPavyiding @NNU by ITNNU and solving for CNNU: 


0.20 


CNNU = .07877 (BI.RE) (ete Gye TNNU. 


ARGPFAC 


Meas form is readily digestible by the computer. #hen CNNU 
Mes been computed the caiculation of one duta point is 


Somplete. 


igen the calculations have been completed for all data points 
the processing of experimental data is complete and analysis 


eer COMMENCE « 








ce 
PNA eee “CONVECTION HeAt TRANSFER RESULTS 


Figures 1 thru 4 are composite plots showing the 
experimental variation of Nusselt Number with Reynolas Number 
for each of the four geometries which have been tested. The 
meriation in slope of the plots as the angle of inclination 
mg@enqes indicates thet the inclination factor By LS Note 


mole multiplying constant but iS an exponential function of 


Reynolds Number. 


Po is defined as the ratio of the experimentally determined 
Peeseit Number to the Theoretical Nusselt Number, calculated 
fSvemere mocitied Grimisonm equation, at a particular Value of 
Reynolds Number. Thus for every experimental point a value 


foe aes Reynolas Number may be plotted. 


le 
Peener than plot these valves by hand, the deta processiae 
facilities at David Taylor Model Basin were utilized to enhance 
miemopecd and accuracy of the process. An existing program 

Wes Modified and used to provide the needed output. This 
ereodram is described in David Taylor Model Basin Report No. 


2037/17} and in Appendix G. 


The program takes as input the values Py and the corres- 
ponding Reynolds Number plus identifying codes which serve to 
individualize each of the points. The output is the equation 
ef three lines of regression which are produced by means of 
meme method of least Squares. In addition, there 1s a qrapiie 
output which plots the three lines of regression which are 
determined. Samples of these output forms are included in 


EOpEndIx Goa complerce sStmmary of the reqression analysis 


coefficients is shown on page 42. 
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Of the three regression equations for each geometry, one 
mepresents consideration of all the points, one is determined 
by all the points within two standard errors, anda the third 
represents the line determined by points within one standard 


error, where a standard error for a set of N points is given 





by the relation: 


2 (v- vave) 


rege | 


= a fader ef AD) POSTE OS ~ 


5S... bans 


ox V 
; where Y avg. is me NI 


Us| 
y 


The choice of which equation to utilize was made through 
inspection of the SC4020 plots of the regression equations. 
In almost all cases, the values of By at low Reynolus Numbers 
were outside the two standard error criteria, and, Since the 
first three points are among those used to make ae first deter- 
mination of exponent, these points influence the phase of the 
regression line unduly. Further, the data for values of Reynolds 
Number less than 2000 is suspected to be influenced by laminar 
flow. Thus only that data which fits within the two standard 


error criteria was allowed to determine the regression equation 


and the seccnd set of coefficients were used for further analysis. 


For design work, the equations as shown in Table I may be 
used 1f the designer is working at a specific angle of inclination 
and geometry. Entering with Reynolds Number as an argument, 

a value of By will result which can be used as an additional 


factor in the Modified Grimison equation. 








For greater convenience, two aesign plots have been 
Meeperea. Figures 9 and 10 in the cnclusion section, give 
the factor By as it varies with angle at tnrec constant values 


Meare ynolds Number. There are l2 curves, three for each of 


me Lour tube arrangement geometries analyzed. 


PReering with the angle of inclination and the geometry, 


the designer may read a value of Py for each of the two values 


of Reynolds Number which bracket the argument he requires. 


By Linear interpol. tion, By Can (be pakeLCularized to ‘the 


desired Reynolds Number. 
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pummeiey Of Inclination Pactor Reqression Coefficients 
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Ror POlLiNcS” Whenin two Stangere Errors 


M 
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0.5087726 LAOS e aL seycu 
0.6856272 4.110530(10)7* 

-0.6487533 0.6318375 
0.8435844 1.62334(10)7% 
0.62085 38 OSS Staelin: 
0.7091857 3.031199(10)7> 


0.7860375 A.310665(10)7> 


1.155682 ~A.277594(10)7° 


7.290824(10)7* 4.853093(10)7- 


3.172470(10)7> 


0.7260181 

0.7262158 Wee any el 0) ie 
0.7653439 1 ODS Sa) = 
0.4970187 Jem ea) = 
0.6867613 Pee Som 1G) = 
Ces0nlTes2 AGO me (io ae 
Oa 751026 DTS Kou ACM OD) ae 


M 
0.6510998 
0.9889 
O.1069 


QO./7754198 


By Shere, 
0.9889 
G2 7eo9 


ie OOD? 


O.2980692 
0.9859 
0.8012598 


0.4712899 


Oe reieys. 
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Mave oLeaOmenn te LONG ouR bot S TANCE 


Bugqueesss! through 66, tnceluded in Appendix H, are plots 
Meetriction factor against Reynolds Number for individual 
Meevetcries and inclinations. Several points pertinent to 
Buosequent analysis are observable from these curves. The 
non-linear behavior of the full logarithmic plot at low Reynolds 
Numbers iS apparent in these plots. This behavior is a 
characteristic of small tube apparatuses that is consistent with 
the observations of previcus investigators and is probably 
caused by laminar flow in the tube bank. In order to keep the 
analysis as general as possible data points with Reynolds 
Numbers less than 2500 have been omitted. Additionally it can 
be seen from Figure 55 that the data for the thirty degree 
meee, S,, = S, = 1.9D 


ab L T? 


all the remaining data and is therefore omitted from the 


In-Line, is completely inconsistent with 


analysis. The inconsistenty is undoubtedly due to faulty 
Bnistrumentation or operator error during the collection of 


this data. 


mie friction factor can be expressed analytically as a 
function of Reynolds Number by fitting the data to an equation 
ee the form: 


PRLEAG eeCONs? (CnRE) 


using the method of least squares. 

A computer program, aescribed in Appenuix H, has been 
utilized to compute these parameters which describe the friction 
factor as a function of Reynolus Number for a particular geometry 
and angle of inclination. The equations which include these 


parameters are summarized in Table II. 
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ie oS ek Od 


0.445318 
0.045465 


0.65543 


UE2ZOLSs0 
0.037337 


O.100758 


0.431623 
Oe S77 oy 
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0.040784 


Uses feoog 
0.074656 
0.017001 


OR OE Srey: 


EXPO 
=-0,457575 


=e GO590 


-0.043053 


-0.081586 
-0.119514 


-0.024007 


-0.0283474 
-0.089713 


=O Ooo 7 9S 
-0.05331 
Os200L2 2 
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Pamemres 5 thru 8 are composite plots showing the variaticg 
Mmmetraction factor with Reynolds number at different flow 
Mages for particular geometries. The plots are based on the 
Meeations listed above. From these plots it may be seen that 


the set for on S, = S, = 1.5D,,, STAGGERED, (2G. 9) 91s noe 


L 
Consistent with the remaining data. The data for this geometry 
was taken by previous investigators who experienced some 


Meeeracuity in correlating it. This previous data will be 


omitted from further analysis. 


mean tie equations BPepresenting Ene VvVarletion of \iriction 
Mecor With Reynolds number, it can be seen that there is no 
Seetectnc pattern in the variction of Gi ther tne coer ficient on 


exponent with angle of inclination for particular geometries. 


Friction factors for Reynolds Numbers of three thousand 
through eight thousand, in increments of one thousand, were 
Bemepuced aS part of the analysis. Program Values of Reynolds 
Number = 3000, 5000 and 8000 are tabulated below to indicate 


the variation with flow angle for particular geometries. 


REYNOLDS NUMBER = 3000 


Pr Leto 4) ae cor 
GEOMETRY 30 A5 60 a5 90 


Se Ss, = 1.5D,,, STG re 0462021 O40 50S 7 Oe 1 140 
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REYNCLDSANUMBER == s000 
mes, =2-5D.,,STG = oe ae 
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Ss =1.5D.; Sp =e~0D 
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Peeeresses friction factor at any angle as a function of tne 


Berection factor at ninety degrees by simply taking ratios of 


Peesvalucs and plotting them: 


where £8 


£90 





Aes 
Ro = £6 


£99 


Ene TrCCELOM LoGlor sch sel algae. 


the £rictiom factor at ninety Gdeqrees as expresseaq me 


EXPO 


£, = CONST (CNRE) 


particuldérized to the geometry. 


CONST and EXPO being 
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The values so derived are tabulated below: 


Nop = 3000 FQ 
SEOMETRY 30 45 60 90 


i, 
Sp=l-5Dep Fs. SG eo Ome oO Je OP aie) shs! 1.000 
Spal mew S,=26 OD,, oiler. Ceo IU0 O.450s 2001S 1.000 


eo 
Nop = 5000 
Sm= Su 1. 2D,, I.L. - OPAL S Haas 2 On 7904 i OO 


Nee = | eh00ks 


Spry = i 2D rent laie aad 0.0455 0.8688 iO 
Spal.oD, S,=2-0D, SIG 0.3470 0.6100 Cla ES Oe i000 


Sp=l.sD., Sy = 2-0Dq, telewO. SLGS2 TO,oz4t 0.6096 JEP AG LOLS 


Values of the ratio FO are plotted and are included in the 
Memerision section aS Figures lil through 13. By entering the 
applicable curve with values of flow angle tie ratio of 
friction factor at the angle © to that at ninety degrees may 
De obtained. For values of Reynolds Numbers other than those 


eabulated, linear interpolation may be used. 
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CONCLUSION 


Convection Heat Transfer 


Figures 9 and 10 represent the plotted-.results of the 
equations developed by the regression analysis of individual 
M@neclination factors. These plots show the variation of 
inclination factor with angle for in-line and staggered 
geometries. Six curves are plotted in each case representing 
constant values for Reynolds Number of 2000, 5000 and 8000 and 
Bearticular values of transverse and longitudinal tube pitch. 
The curves can be used directly as a design aid by applying 
the value of inclination factor obtained as a multiplier to 
the modified Grimison equation. Values of the inclination 
factor for Reynolds Numbers other than those plotted may be 
obtained by linear interpolation. The modified Grimison equation 


now becomes: 





e 


iO. Coe ree ee En (Np : 


0.260 
DA oe UD ( 


oa 


ie Ssnould be noted that the variation Gt anclinetion factem 
with angle of inclination is by no means linear so that there 
are regions for both staggered and in-line geometrics in which 
relatively large changes in angle of inclinution cause small 
changes in the rate of heat transfer. The reverse situation 
also can be obtained. The choice of which region is optimum 
must be left to the user who is familiar with the requirements 


peoetimitations of his particular job. 
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Pere llONAL RESISTANCE 


ierojurecm el through Ms represent tne plotced results 
Mamptne analysis of frictional resistance for each of the 
Myree geometries for which reliable data was obtained. Each 
mmeere Contains three curves, cne cach for values of Reynolds 
Hamoers cof 3000, 5000 and 8000. Additionally on each figure 


Mere appears an equaticn in the form: 


f = (Fa) woe 


ines 1S the equation for the friction facter in each casem 


Two uses of the curve are suggested. The first is to 
@epculate an absolute value of friction factor. This method 
should be used with caution for while the method of obtaining, 
reducing and analyzing frictional resistance data is believed 
to be more accurate that that used by previous investigators 
it does not agree closely with their work. The second use 
a5 to estimate qualitatively the effect of varying the flow 
meee On the frictional resistance in a particular piece of 


heat transfer equipment. 
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PM COMMENDATIONS FOR FURTHER INVESTIGATIONS 


lemme poeraktus, teSteprocedure and method of reducing 





and analyzing data in this thesis represent a suitable method 
of empirically determining the effect of flow angle on 

@envection heat transfer and frictional resistance in a tube 
bank. The results presented are limited in that they include 


Bmery four possible tube geometries. 


It is recommended that the series be extendeda to incluae 
all values of transverse and longitudinal tube pitches likely 
to be found in modern heat exchangers. Such work would involve 
manufacture of new tube banks and sealing strips. It is 
PUrth. r recommended that data pe taken at an angle of seventy-— 
five degrees both for other geometri:.s which may be tested and 
for those for which this angle has not been previously used. 
The additional point wcula be ot great assistance in firmly 


Meecting the upper end of the inclination factor curves. 


If a blower of larger capacity should become available 
the test data for all geometries should be extended to higher 
feemelas numbers. Prospective investigators ere advised that 
eee limit of the currently installed pressure instrumentation 
Nas been reached by the present air supply system. They 
ere also aavised that the manufacturing processes are long 
Seas cedicus and that the eccondition of the Laboratory steam 
Plant should be carefully investigated before unuertaking to 


use lve 
che reduction anavanalysis of test data cculd be facili Glee 
by combining the three cumputer programs usea in this work into 


Be progrem to reduce and analyze data in one computation. 
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APPRENDIA A 


Appar acus 


At the onset of tne project, it was felt that the apparatus 


as designed by Bond ana Wallin and later used by Kiss and Smith 


would not be adaquate for tne task without considerable 


Meaitication. The apparatus as constructed and available had 


been used for one basic geometry by Kiss ana Smith and had 


an inherent disadvantage in that it was difficult to make any 


alteration in tube geometry. 


gme Gesign of a new test section was therefore undertaken; 


The principles which guided the design are: 


1. 


The section must be of the same cross sectional 
dimensions as the original apparatus aS verified by 


Bond and Wallin. 


The section must have the capability of rapid tuce 


geometry changes. 


All Manufacturing eperatLiens mecessery for completicg 
of the design must be feasible with equipment avail- 


able in the Webb shop facility. 


The design must utilize the same tube banks that were 


prepered for the previous investigators. 


A design which satisfied the four st.ted requirements was 


developed. The configuration is generally the seme as that 


used in the earlier experiments. The steam supply and drain 


line were used without moaification, as was the air supply 
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meeeemn including the orifice meter. The entire test section 
was rebuilt in accordance with the plan shown as Fig. 15 in 


this appendix. 


Mine eest Section Consists of two rectangular wooden boxes 
closed at one end, which serve as the entrance and exhaust 
mertei1ons of the test unit. Connecting these two end pieces to 
morm 2a Continuous square duct, is a center section of aluminum 
alloy. This center section is designed with stationary side 
plates and movable top and bottom plates. These features 
enables the tube bank tc be adjusted to any angle from verticel 
through 60 degrees from the vertical. This was the practical 
Beeent of adjustment in that the tube banks were not long 
enough between headers to permit another range. This is, 
however, considered to be a sufficiently wide range, as it 
permits angles of flow incidence from 30 degrees to 90 degrees, 
at any practical increment. This design parameter sized the 


mength of the aluminum center section at 30 inches. 


The end sections are made up of 3/4 inch exterior grade 
Marine plywood. The joints are made up with glue and joined 
Pears x 1% " wood screws. The sides of Ehe section are 
routed 3/16" by 3" to receive the sides of the aluminum center 
PeGceron such that the center plate wili £1 flush. The 
wooden sides are also routed at the top and bottom edye, 3/16" 
wide by 4%" deep to receive the edges of the sliding top plates. 


The top and bottom of the wooden intake and exhaust sections 
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Meeerouted to a depth of 3/16 inch to effect a flush fit for 
the sliding plates of the aluminum center section. They are 
also routed across the width and to a length which will 
accommodate the sliding plates through the range of incliné&tion 
Memeene tube bank. This range is shown in Figure 17, the 
Meawing of the section. The intake section as a 6 inch 90 
elbow secured to the end for delivery of air through a 6" 
diameter hole. The exhaust section has a similar stove pipe 
elbow made up to the side of the section. The heated air 
mmmpeausts through a 6 inch diameter hole in the box then through 


the stove pipe exhaust line to the atmosphere. 


Both intake and exhaust sections were sanded smooth and 
finished with two coats of shellac prior to assembly. After 
assembly, all edges were caulked with a silicone rubber sealant 
and the stove pipe connections were made up with asbestos cloth 


gasketing using Permatex compound for a leakproof joint. 


The center portion of the test section is fabricated from 
2024-T3 aluminum alloy. This alloy was chosen for its free 
machining characteristics, since several machining operations 
were required in manufacturing the numerous sealing strips 
between the tube banks. It was later learned that the temper 
of this plate was not ccrrect in that when the plates were cut 
Co size and the edges and sliding grooves milled, the plates 
foamed OUL Of planevduc ce locked Im rolling Stresses. Arter 
several unsuccessful attempts to straighten the plates by 


different methods, the decision was maue to anneal the plates 
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maenolates and thus streesrelieve them. Several industrial 
contacts were made, ana threugh the courtesy and kindness 
Mmmpene Sylvania—-Corning plant in Hicksville, N.Y., the plates 
Meee annealed. Proviaged thus with plane surfaces, construction 
of the center secticn may proceed. One further note could be 
added for the benefit of those who may be curious, a more 
Mmeoper alloy choice for this epplication is 2024-T53 this 
temper being stress relieved by stretching. The side plates 
were made up with all edges milled and fastened to the wooden 
end pieces by flat heat wood screws, courttersunk to effect 

@ plane surface. The top movable plutes were mauve up having 
all edges milled and two 3/16" wide by 1/16 inch deep groove 
milled in the sidés to effect a seal with tne sige plates and 
—O provide a channel for guidance of the moving plates. Each 
top and bottom plate (four total) has a projection along each 


edge which slides in ae channel. 


The channel was mage from an aluminum storm window shape 
fect oved to have insufficient strength to hola the weight of 
the tube bank and secure the plates in position as well. To 
correct this, heavy steel yokes were made from angle and 
threaded to close the section and secure the sealing strips in 
place. All but one of the channels were removed, the one 
Bemalining only to serve as a positioning device. Manufacture 
of the sealing strips represented the largest expenditure of 
eime for any one operation. For each of the 12 combinations 
of angle and geometry, 22 sealing strips were required, lil 


fech tor top and bottom. Each sealing strip required edge 


_ 














67 


fmeing, groove milling, and 22 holes arilled or milled to 
close tclerances for compl: tion. In all, 5680 holes were 
machined, and a total of 264 strips were preparea. The sealing 
Strips proved very satisfactory in maintaining the tube geometry 
while providing excellent sealing characteristics. The only 
modification required was to provide a small hole in the end 

of each strip to permit the separate pieces to be wired to 

hola them in place while changing tube geometry. An isometric 
Meeechn Showing the arrangement of the center test section 

and the sealing strips is provided as Figure 15. Drawings 
Piiewing the sealing strip as a part of the test section and in 
M@etatil appear as Figures 17 and 18. A photograph of the 


Sediing strips in uvetail is available as Figure 22. 


Baer seccion 15 supported = tromna laboratory veable by twe 
wooden frames and a framework of steel angle which was modified 
ey tO facilitate the changing of the geometry. This 
modification consisted of welding an angle on each corner of the 
metal structure to proviue support for tne tube bank headers 


when changing the sealing strips. 


The tube bank is thet used by previous investigetors. 
Tt is made up of 10 rows of 4% inch O.D. copper tubes, 22 tubes 
her row, with each row of tubes being fitted with an inlet 
and exhaust header of 14 inch ccpper pipe. A detail drawing 
of an individual tube bank is shown in Pigure 16. The inlet 
headers are connected to the steam supply moisture trap by 


3/8 inch copper tubing. Many of these lengths of tubing had 


_ 
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had to be replaced to permit the freedom of movement required 
in the tubve bank. The exhaust headers were connected in the 
same manner to a drain header of 14 inch copper pipe. An 
adagitional problem was that these lines kinked when movea about 
due to strain hardening and Several were replaced for this 
Mmeason. For some geometry changes, the exnaust and supply 


lines were completely disconnected to avoid this problem. 


Two small holes are provided in one side of the center 
Mest section for thermistor leads. Four holes are provided 
Merepressure taps in the test section. Two are in the inlet 
SBection, positioned in the middle of the vertical height and 
menenes from the center section. Two taps are provided in 
the exhaust section on a similar vertical positicn, 4 inches 
from the test section. Thermistor probes are provided in the 
6 inch inlet piping just before the entrance elbow and in the 
Peinel Outlet section to avoid radiation effects from the 


tube bank. 


Previous investigators had used various methods of sealing 
mee bonks; Keferences 4 and @. The method used in this series 
Peetests, though long and painstaking in preperation, proved 
Worth the effort while testing. There was no need to use any 
Sealing compound cnce the tubes were in place and the top 
and bottom moving plates mace tight. The tolerance to which 
the sealing strip holes were prepared was 0.008 inch and there 
was little or no leakage. Further, a geometry change could 


be made in about two hours on the average and there was no glue, 











isecaier, etc. to clean from the tube to cause further delay. 
In this regard, the extra effort was rewarded and this 


arrangement is recommended for later investigators. 


Figure 14 shows an overall isometric view of the steam 


meant end test apparatus. 
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Overall View of Test Section = 90 Degree Configuration 








Fig. 22 Sealing Strips. Configuration Shown is Sy= 2.0 D4, 
Sts 1.5 De, Staggered Tubes, 60 Degree Inclination. 
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Tnstoeumentlation 
Hctipemoacure 


Instrumentation for obtaining temperature readings was by 
Meermistors installed in air intake and exhaust lines ana on 


mubes to sample temperatures throughout the tune bank. 


Thermistors are solia state electronic semi-ccnductors 
Miach exnibit a large negative temperature coefficient of 
resistance. Their small physical size and the large negative 
Me@oereture coefficient of resistance make thermistcrs 
perticularly well suited for a, plication in small scale exper- 
@m@ental models as units can be locdted within a very limited 
Space and great accuracy can be obtained with a relatively 


Semole electrical circuit. 


The thermistors used were Fenwal GB35J1 bead type thermistors. 
Mhese units are 0.043 inches in diameter and have a nominal 
Beetstance of SK ohms + 20% at De Game  tMermisters were 
Calibrated by the use of a constant temperature bath. Bolling 
Water and boiling benzene were used to obtain calibration 
Deyonad tne ranged the bath. A sample calibration curve 1s 
Maeliuded as Figure 23 at tne end of this section. Thermistors 
MWsed for inlet and outlet temperatures were fitted into 
Mmeeoes Made from brass tubing and sealed with silicone rubber 
Meeeeant,. ihe probes were located in the center eerene inlet 
Meee xnaust ducts atl a point around a nanety degree bend from 


the tube bank so that they would be screened from direct 











radictiom from the tube bank. Some difficulty was experienced 
my atifixing the thermistor units tightly to the ture bank se 
Mie a true tube temperature would be inaicated. The problem 
was finally overcome by taping the thermistors to the tube 

fee plastic electricians tape. Prior to placing the 

thermistor on the tube, a small groove was ground on the tube 
Surface to receive the beaad. The tape and pigtail leaas were 
then covered with a high density epoxy. This effectively 

held the therm:istor in place even during tne curing period. 
Thermistor surface sensor assemblies consisting of thermistcr 
Beaas bonded to teflon or silicone impregnated fibergless tape 
are commercially availavle. Tne cost of these sensor assemblies 
iS approximately eignt tim s that of tne thermistor units 

so thet it was not considered economically feasible to use then. 
FOr Similar reasons it was decided to undertake calibration 
Beethe units rather than to purchase units which track a 


meaenGard Curve wiLEhin nerrow limits. 


The thermistor leads were connected to a selector switch 
which switched individual units into a wheatstone Bridge circuie@ 
The resistance necessary to balance the bridge was measured 
ena this value was used to enter ’the Gelibration curves ve 
determine temperatures in aegrees Farenheit. A diagrammatic 
sketch Of the Cempereature measuring Cireiul 2 appears as 


Figure 24. 


A high failure rate was experienced With the chermistce 


Pits. Ine ha qnest perGcemlage Of temiiine. Wereuaic sua 











er 


accidental severing of the small diameter wire pigtails at the 
point where the wire entered the glass beaaq. These failures 

were caused simply by the rough service to which the units 

were subjected, particularly those installed on the tube 

surface. In addition, several thermistors were found to be 
mmegcic 1n behavior, culminating in a short circuit within 

the bead after a short time. No explanation for this performance 
has been found. Faulty units were replaced with spares until 

the supply was exhausted. At this time, it became necessary 
Mmemcreplace the exhaust temperature thermistor with a mercury] 


in-glass thermometer. 


As an aSide, it should be added that the thermistors were 
Mee Given fair trial in this application. Time necessity te 
move the tube bank about roughly placed undue strain on these 
Mme s. based On this experience, hardy commercially prepared 


mounts such as those mentioned above are recommended. 











Pressure 


instwgumentation for pressure was by manometers containing 


Oll or water aS appropriate. 


The manometer used to determine the pressure on the 
upstream side of the orifice was an open water manometer. The 
barometer reading was converted from inches of mercury to 
inches of water and added to the manometer reading to obtain 


an absolute pressure upstream of the orifice. 


An oil filled U-tube differential manometer was used to 
indicate the pressure drop across the orifice plate. A water 
filled U-tube differential manometer and a slant tube manometer 


were used to measure the pressure drop across the tube bank. 


A Bourdon tube pressure gauge and a bellows type compound 
Gauge were used to meaSure steam inlet and drain pressures 
respectively. 


Flow 


Instrumentation for weight rate of flow was by thin plate 
orifice meter with pressure tops located at one pipe diameter 
preceeding and one half pipe diameter following the inlet face 
of the orifice plates. To obtain@a large spectrum of Reynolds 
Numbers, three separate orifices of 2 inch, 3 inch, and 4 
inch diameters were used. The weight rate of flow was obtained 
from the differential head across the orifice. It was found 
that it was possible to cover the entire range of weight rates 
of flow which could be obtained with the available steam 


~urbo=—blower with the 2 anen orifice, 
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Operation of the test equipment can be divided into two 


Metegories: 1) Lining up, lighting off and operation of the 


steam plant and auxiliury equipment and 2) Operation of the 


mest apparatus itself. 


The following procedure is followed in putting the steam 


plant into operation: 


1) 


2) 


So 


4) 


a) 


6 ) 





The cooling tower is checked and filled as necessary 

BO CYSULe Ha “OrOper GUdnt. by Cl eocling water, 

Cooling seeveo met ais vd edema l. CUla rite suing care 

Seer bea. 

The “eemdensete Pulipewsestamred. 

Fires are lighted in the boiler fcllowing the procedure 
given in the boiler cperetion and maintcnance manual. 
When the steam pressure has reached abcut 100 psi the 
air ejectors are lined up and sufficient nozzles cue 

im CO mMaintall obour ten anches of VecuuMoin tie 
condenser. The most satisfactory operation was obtained 
when both second stage nozzles and no first stage 
nozzles were used. 

The turbo blower is lighted off using a throttle pressure 
of about 5S psi. The blower lubricating oi! tank shoud 
be checked to ensure that an adequate supply of 


I blelanttes=enkie(s) out vals) salle Saiaverr 








7) The steam separator is drained ana the tube bank 


steam supply and drain valves are openea. A pressure 
Of ten psi dat the steam Separator was found to be 


VELyY Secisiactory. 


When the above operations have been completed ana the 
Menaitions in the apparatus as reflected by inlet, exhaust, 
and tube temperatures have stabilized the procedure enters its 


second phase, the collection of data. 


Mest procedure consists primarily of varying the quantity 
Se amounts of air supplied to the apparatus and recording the 
Mera provided for by the installed instrumentation, as dqescribed 


in Appendix A. 


The weight rate of flow of air to the test section is varied 
eecirottling the blower supply steam with fine control being 
Pecomplished by manipulation of the aqischarge vleed valve. 
Poacisfactory operation is obtained by starting at low flow rates 
mide Lncreasing graaually to high rates. with this precedure,; 
mlereabparatus stabilizes quickly so that the amterval setween 
runs is very short. It should be notea that because the test 
section and the turbo blower are supplied by steam from the 
some line fluctuations in boiler pressure caused by variatione 
in blower demand will cause fluctuations in the pressure to the 
tube banks. As the tube temperatures and exhaust temperatures 
ere Sensitive to the supply pressure Variations should be 
minimized by ensuring the tube bank supply is the same before 


each run. 


4a 


i 


_ 





The pressure droov across the orifice is indicative of 


[row rate, accordingly successive runs may be maue by adding 


small increments of pressure drop by throttling the blower. 


DURING Vees ling, “lest somaevsneageous to Make a“runring 


mrot Of the quantity TT-TI Y DHF versus\/ DHF on full 
Bid WE 


iegarithmic paper. Date / DHF may be considered 

TE 
mirectly proportional to Nusselt Number while\/ DHF may be 
considered directly proportional to Reynolds Number. A linear 
plot on full logarithmic paper then serves as a check on the 
consistency of the data. Care should be taken not to take 
figem a plot aS indicative of the valiaity or accuracy of the 


Meca as the proportionality is not exact. A sample plot appeam 


as Figure 25 at the end of this Appendix. 


Orifice plates of threc dififerent diameters are availlap me 
Mer use with the asparatus. Utilization of the three inch 
Orifice permits data to be taken over the entire range of 
Ptow Obtainable with tiie turbo blower. The othér two orifice 


Sizes were not used. 


If an immediate change of geometry is anticipatea upon 
Sempletion of data taking, allowing the blower to run ateer 
the steam supply to the tube bank has been secured cools the 


metal portions of the apparatus quickly and facilitetes hand Tage 


Securing procedure is essentially tne reverse of the 
Pigiting off procedure. @eUpom, securing im celia weather the 


cooling tower should be drained to prevent freezing of the pipes. 
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AEP ENOEX Cc 


Development of Equations Used in Computations 


Computerization of the calculation procedure required the 
development of equations for a number of quantities which are 
Msually tabulated or taken from Curves in engineering 
applications. Development of equations which are applicable 
in the range of temperatures and pressures which exist in the 
apparatus reduce the complexity of the program and permit a 
great saving in core storage within the machine. The quantities 
for which equations are developed herein are: 

toe Thermal cCenduceivity of air 

Py Dynamic vascosity Of air 

3) Specific heat of wet air 

4) Flow coefficients for orifice meters 

5) Expansion factors for orifice meters 

6) Grimison Arrangement factors 


7) Nusselt Number by Colburn equation 


The symbols used in this appendix are those developed for 


use in the computer program. 


1) The thermal conductivity of air varies in an essentially 
linear manner in the range of temperatures which occur in the 


apparatus. An equation may be written in the form: 


k = (k 
nT xf = eC 
Vvhere 
K = thermal “ECcnduceivity of air’ in Beu 


hereto f 


Ak = Change VWniseehermarmecnouctivity per eqree ar cnileimis 
eek 








ay 


Temperature in degrecs Farenheit 


S a constant 


tl 


Values of thermal ccnductivity used to develop the equation were 
Meeon trom Table 2 of Gas Tables by Keenan ana Keys. [9] Tie 
Meation for thermal conductivity is: 
k = 000024 (T-40) + .0144 
Og 
k = U.000024T + 0.01344 
This equation gives values which correspond exactly with the 


tabulated values between 40 and 240 degrees Farenneit. 


@) Dynamic visecsity, like thermal cenauctivity, vVories 
linearly with temperature in the range of application and can be 
expressed in the same form: 


: 2 34 SC 
i 
“where: 
pe = dy ndMase Vi sCosi cy ia Los 
sec-ft 
EEG e = e . . > e 
Se = CN aNGe fi Oy cin eave Gos Ike, Noe rma C meeme ca mOeM ie i.e 
aT 
i = lenperatume ttieaG@egeees o) omenuweie 
: C=) Cons lance 


' Meiues used in the Gevelopment of tie equation were Cakenm feeq 
Table 2 of Gas Tabl.s by Keenan and Keyes. The equation for 


Synamlc Vi —cOSity 1s; 
a 


| 


| +26 (T-40) + LL: | eee 
or 


E= F 112 | 6 


AA 


-7 


Jt 

















oa 


3) The development of equations for specific heat of 
Pe— air @e constant pressure requires the selection and adoptaam@ 
oe applicable equations for the specific Neats of air end water 


and combining them with proper weigi:t.ng. 


The equations selected were taken from Table I of 
Thermodynamics by ecass [vl 


mor dry alr: 


eon) = 021 = Olsae meer = oon x Noe 
and for water: 
een = 1102 See oe 416 
ap 
where: T = average temperature in the tube bank in degrees 


Rankine the specific heat of wet air is then determined by: 


CPw = SPHU (CPH) + 7000 CPD 
Sonus 7oCe 


where: CPW = Specific heat of wet air 
GP = Specific Neaepotneny al: 
CPH@= Speclile ww eareoumwarc™ 
SPHUG=e Speci lie unworn oliow eer pound 


4) Flow Coefficients 

For square edged thin plate orifice flow meters with taps 
one pipe diameter upstream and one half pipe diameter downstream 
Ene following equation for flow Cocif tcientiwovrains. 

Peas Ont ae 
Wnere: 

AK = i low CoecEtieG@ent COrresponatnG foal y ypartrecu var 
Sel Or the parameters pipe Gdliameter KD) ga reatlo vot ‘ratice 


diameter to pipe diameter £8,, and pipe keynolas number (ANR:E) 








Ko = Flow coefficient when pipe Reynolds number becomes 


infinite 





b@eran empirical multiplier 
ene C0) 


V ANRE 


Ko and b may be expressed in the following forms: 





imo =) (0.66014. — 0.01252 p74 Sa) Bowie £0 725 / pv 
0.00025 
p°f* + .oo25p + B* 41.5418 
b= 


(O.0002 + 0.0011) + (0.0028 3 80.0004) 
D D 


|B° + (16.5 + spar) 


POL slCudtilons "SUCh eas Cx2 St. Wn wee veesctE vorpal aculs. Wiese 
there are a limited number of possible values of pipe diameter 
and diameter ratios specific values of Ko and b may be 
@aiculated and applied to the basic equation. It should be 
noted that although 91° represents a number Less Chan unity 
raised to a relatively large power it cannot be ignored if 


meceptable accurecy 12S to be obtained. 


Mae values fer Ko ana b fOreehe Varlouc Oorarice pieces 
avallaple are: 
TWO seh Jor iG weer 


Ko Osco 7 oe 


I! 


D ~<COCED 


Enree Poche, mee. 


ko Os ol SG 


tN 


b O20 ise 


POUrl Ench Cor erace: 


Ko C5075) 


._ b 


O200229 








aes) 





The expressions for flow coefficients for each orifice then become: 
Twominch Orifice: 


hKw=O,59 7310.00 


——— 


ANRE 
three inch OfFlrTice: 


AK = 0.6180 a eS 


eee eee 


JANRE 
Pour inch Orlilce: 


AK = 0.6735 + LED 
JANRE 


Ss Expansion Factor 
Tne equation for the expansion factor through a Ssqueme 
edged thin plate orifice meter witi taps one diameter upstream 
and one half pipe diameter downstream may be expressed as: 


Viocll = (0.4140, Saree lee 
K 


where: ear tact lO Of sOriticerd 1 ane tem sCOupl pe wa iane ter 
XA = ratio of differential pressure to absolute statie 
inlet pressure 
K l=" atic Of specie ve Meate. 
Lo a COonSteantteVelUGROt gy) te oe SSC ian ey euuice of Jam 


ere calculated for the two, three and four ztnch orifice meters 
the equation reduces to: 


FOr 2c OrlLticess YY =e) aS pe 9,8 


FOr 23) Orit ice te aoa S30 7x 


For 4” orifice ey e-— le 520% 
POrCULIIIZaLIGn an tne ~Logram the ectorm. muse be 
expressed in the same Units as the input data. The numerator 


X 1S equal to the differential pressure across the orifice 


4 


4 
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meter (DHF) which is in inches of oil. The denominator is 
moual to the sum of the barometric pressure (PB) expressed 
in inches of mercury and the pressure above barometric ahead 
of the orifice meter (P,) expressed in inches cf water. All 
units will be converted to inches cf water so the expression 
moor xX becomes: 

P= Creoe Orn 

Py spel Shs spose sls: 

where: DHF, Py and PB are i1nSue Gqudmerties as defined in the 


preceeding paragraph. 


0.827 = the specific gravity of the oil used in the manomewes 


ie S0G= Lie COnversiomet actor sbom 1 mches Gr MeLreury 


Co. 1nChes- GEewatelr. 


Inserting the expression for X and simplifying the expression 
Ber expansion factor as used in the celculation are obtained: 
MOR (aot id Ce . 


ee a 0.24464 DHF 


Py rie le Pye he] 9 Sais} 


mor 2o" OFrlvilce: 
Yeoit)f]le- OO. 2525 7/7 Doe 


Py + 13.596PB 


Mor 4” ofitiec: 
Yedl- Seis sb ague 
Py oul oi 6) Sepals. 
6) Grimison Arrangement Factors 
The arrangement factors used in the Grimison Equation 


ere a function of longitudinal and transverse pitcnes, Keynotes 


Number and the arrangement (staggered or in-line) of the tube 





Se 





Pete, ana dre noc amenable to a simple analytic expression. 
Mowever as only one transverse pitch, 1.5 tube diameters, and 
eelimited number of combinetions of back pitch and arrangement 
were tested it was found feasible to plot arrangement factor 
as a function of Reynolds Number only for each geometry. The 
Beeves SO plotted were parabolic in form. An equation of the 
form C, (BNRE)* + C, (BNRE) + C, = ARGFAC could be written for 


each of them by utilizing the LaGrange Method. The equations 


memacveloped are of the form: 





ARGPAC = Cy (BNRE)* ae C5 (BNRB) + C. 
where: 
= LO 5 
GEOMETRY Cy x (10) C5 1 OD C. 
Sm = Sy, an as) D, 
IN-jLINE woo -.2963 JR OISSS: 
STAGGERED ~lec? Bote) 2 94059 
Sin = 2a Diy Sr =e D. 







IN-LINE 
CaO 


WY 
UW 





Om = 
SLAGGERED 


eo 






Nusselt Number by the Colburn Equation 


in order to establish a lower limit wnich would be comoveraga 
me warallel flow cae Colournm equetron was used to calculate 


Nusselt Number. The Colburn equation is: 





me ee GD, (spear 
R == | 
BNere: 
GD, = Reynolds Number based on hydraulic diameter? 
{4 


For the Colburn Nusselt Number to be airectly comparable with 


thcse developed experimentally or by the Grimison equation it 








96 





must be based on tube diameter rather than hyaraulic diameter. 
memaccOomp~ish this both sides of the eaquution are multiplied by 


the ratio of tube ui.meter to hydraulic diameter. 











OR S16) 
eet as Saene 3 en D. cin y hee? 
K D ee a as cas 
H yy a, 
h_-D 0.30 OeZ20 ae 
Se one Ge ae CPh 
k t - mat MA 
ae rT) k 
ry H 
in the ncmenclature used 
' ,,0.80 CPT 6 - Nee 
CNNU = 0:;023 (BNRE) De (spe) 
ae 


Mere nyaraulic diumeter 1s defined as four times the free flow 
area divided by the wetted perimeter. For a staggered tube 


eerangement: 








2 i oes 
a | SpS, ya 7 
| Ws) ‘a }- ae eee 
\ \ ; 
D a Ze 
yw Le J | | 
: 2 x | ; 
\ Di 
: fay Cc = 
a ae = eee e ae | 
Seer 
Benet LEE Ran eee ee 
vi Dy r 


2 
= as: = 
Di = 4 | pr TT On 
i De 
Pe Z 
m2 
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Since only one value of Sip and two values of S, are considered, 


L 
the hydraulic diameter can be numerically evaluated for these 


values. For S, = 1.9D,, S, = 1.5D, 


D. = 4 (1. 5D,.)* ~ ptt 





9! ts 
Be In! 
D., =| 9— TT D 
H ———— te 
(25 : 
eemilarly for St = 1.5D,, Dig= 2.0D, 
aay 
The factor Dy eae then can be evaluated for a particular 
oh 
geometrye Thais factor has been given the name DIAPAC anGdiam. 


mraeluded in the Colburn equation to give it the form: 


CNNU = .023 (BNRE) 2°89 (DIAFAC) (sey 
k 











APPENDIX D 





Experimental Duta Reduction 
Progam Descripeion 


Piece cemellten progr atmaeser Vedmmeroim 25 c PORTRAN LV 
program designed for use with the IBM 7090 data processing system. 
moe program calculates Nusselt Number, Reynolas Number and 
Mmerctional resistance coefficient fdr varying flow angles from 
Beeerimentclly derived data in the Haeberle Laboratory, «ebb 
Institute. The constants ana to some degree the method used are 


Meeticularized to the apparatus used for data collection. 


Mne first step in the program 1s to calculete the sines ves 
the angles of inclination at which tests were conducted ana to 
asSign to them flocting point variable names in order thet they 
may be repeatedly called without recelculation of the sine 


Pumction.e 


Comtrol of the program ds acileved tnrougnea YO loco 
utilizing the variable M which represents the sequential number 
@roa set of gata being processed. it is required that the 
maximum value of M called KASES be sup; lied as input each time 
the program is used. The value of M's incremented by one each 
mane ad Set Of Gata #5 processed. The UC statemene refers cic 
mechine to a READ statement causing a New See Of Gata to be 
read in and processed. when the maximum value of M (M= KAS&£S) 


ms reached the program 15 completed. 
Required input to the program is; 
dye Runi code emumbe: nme Our Cigiteda xed epOlntVoriablce 


(IRCN) 














2) 


3) 


4) 


5 ) 


6 ) 


ae) 


8) 


a) 


Oo 


PA.) 


12) 


ine, 






ao 





Pressure drop across the thin plate orifice meter 

i eenenesS ) OL Odueea Cirecsonu TOUT Gigit floating polms 
variable, (DHF). 

Pressure above atmospheric ahead of the orifice meter 
in inches of water, a three or four digit floating 
pount Vari dbte weur sd). 

Barometric pressure in inches of mercury, a four digas 
floating point variable, (PB). 

Pressure drop across the tube bank in inches of wateze 
a three or fourveigitsflocting point Veriable, (Seye 
Temperature Of ene Inlee vale sin degrees Ff arenes ta 
three digit floating point variable, (TI). 

Temperature of the exhaust air in degrees rarennheit, 
asfour digit floating pointy Varicble- sts. 
Temperature of the tube surfaces in degrees Farenheit, 
@ LOUETGL Gilt Lieltang poem aya uunaiolc oma l Je 

The diameter of the orifice plate in inches, a single 
digit fixed point variable.) (id). 

The-specific humidity of the Wnlet air an grains, pes 
pound, a three digit floating point variable, (SPHU). 
A geometry factor indicating whether the tube geometry 
is staggered or in-line, a Single digit fixed point 
variable, (LGF). 

An angle factor angicating the angle Gf anclinationmmes 
the tube bank, a single digit fixed point variable, 
CLANGLE 7. 

Back pitch factor. indicating wnetmerwsine back pitce 


is 1.5D or 2.0D, Single digit fixed point variable (Ua 








ie eeOGr et seme eULaLeSs Ene Mitermediare quanti Gy 


Seftor us@ein the pipe Reynolus Number calculétion. 


The program then queries the value of orifice diameter to 
@etermine the path to be taken in calculating the Reynolds 
Memeer thru the orifice. 4 Tne Constants chiesen tcr Celculatwume 
ene expansion factor and flow coefficients are functions of the 
diameter of the orifice meter. The procedure followed is 
moe Same for all orifice diameters. A value of flow coefficiene 


(AK) is assumed, and a value of V is calculated. 


The three quantities which have been calculated are 
Mimeciplied together to determine an initial estimate of Reynowam 


Number. 


When an estimate of Reynolas Number has been made a better 


estimate of flow coefficient is obtained. 


The estimate of flow coefficient is then compared to the 
Material value. If the two values vary by cone halfz of one per 
Peiemor LeSse tne Calculeareamalve Of sewed a.) eclce necioemeteiie 
of the flow coefficient and the value of Reynolds Number 1s 
eeeeptead and the calculaeron Continues. i12f “che ditirerence 
metween the two values of flow coefficients 1s greater than one 
malt of one per cent the calculated value of flow coefficient: 
ms used £Or am Enilielewalve Gor, TecelGulocuiommouse eye pac 
Nmumber end) fiew coefficient. This process Coneimucs Until eae 


ene Nalf of one per Centecriterion 15 reached. 


When Satisfactory values of Reynolds Number ana flow 


coefficient have been established the program proceeds by 





cep! 


Piering ene Orifice diameter and Converting it to a floating 
Mellie Variable fOr USe Pmethe weigne rate of flow calculatioug 


Weight rate of flow is then calculated. 


The average temperature in the tube bank is calculated 
and converted to degrees Rankine. With this average temperature 
the specific heats of dry air and water are calculated and 
combined to determine the specitic heat of the air in the 


BODAaractus. 


Log mean temperature difference and film temperature are 


then calculated. 


All quantities necessary for the calculation of Nusselt 


Number are nbw available and line Quan Eley s ea eulacted. 


The angle factor (LAINGLE) is used to select the proper value 
Sr the sine of the angle of inclination previously calculated 
This sine function is then multiplied by the calculated Nusselt 


Number to determine the Nusselt Number for the angled configuration. 
The Reynolus Number in the tube bank is calculated. 


The Program then queries the geometry factor to determine 
if the geometry is staggered or in-line and branches accordingly 
RO calculate film temperature for the frietion £actor 


Sompucation. 


Reynolds Number for use in the friction factor calculatzem 


Pos Gqalculated. 


The density @€ the air in the apparatus and the friction 


factor are calculated. 





ST fey tay 
Lt} 2 
et TE ar Ja te 





After the friction factor has been computed the program 
Poeain queries the geometry factor and then queries the batk 
pitch factor to determine the geometry so that tne proper 
Merangement factor and diameter factor for use in the Grimison 
and Colburn ‘equations can be computed. when the arrangement 
factor h:s been selected anu comcuted it is employed in the 
Grimison equation to calculate a value of Nusselt Number. The 
value of Nusselt Number calculated by the Grimison equation is 
Mrvyidaced into that obtained by use of the experimental wuata. F 


metlue of Nusselt Number is calculated using the Colburn esuationg® 


All computations have been completed so that an output 
ieipe May be written. Output consists of: 
1) Run Code Number as read in (IRCN) 
2) Nusselt Number determined from experimental cata (BNNU) 
3) Nusselt Number determined from the Grimison equation (T N 
4) Nusselt Number determined from the Colburn equation (CNNU 
5) The ratio between the two Nusselt Numbers (FTHETA 
6) Reynolds Number based on the film temperature used in 
the Nusselt Number calculation (BNRE) 
7) FriCctgonefractor ma REP AG) 
8) Reynolds Number using the film temperature used in the 


fret One Geter see Leu oe tom, 


Upon completion of writing the outpuL data another Seu am 
data is read in anu the program repeated. With tne 7090 system 
calculétion of each set of aata _equires about one seventh of 


one second. 











yh Dleek GlaGgram Of 9Ewewproqram 1s presented herewith. Fae 


Merock diagram is a graphic representation of the data processing 
system procedures used in the solution of the problem. The 
diagram is inciuded to assist the reader with a visual record 
Of Operational and logic steps in the program and the relationship 
meen step in the program has to the solution. 

The symbols used in the diagram are the standaid IBM symboue 
representing machine functions. A Key to these symbols 15 
agncluded in the Nomenclature section. 


Pe POGran cuisting alsent ol lov saene blog Gicgran. 


The computer facilities used for these computations were 
mmesec at the David Taylor Mouel Basin. Data was prepared using 


Paguipment available at Webb Institute. 
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Experimental Data 


cil Ss we penGieemereCoCi@s, Cums tingwonralim test data, line 


data is listed sequentially by run code number:. No lcbels 


meemeLovaided on the data listings, accordingly the following 


foacormation is provided: 


Solumn 


if 


1Q 


Jal 


12 





DeSénipcion 
Run code number for data point identification, 
Pou giao nem a 
Pressure drop across the thin plate orifices, 
incnes ger. Ol | wen reewOl (hOUEeO1gits {DHE ) 
Pressure ahead of the thin plete orifice, imei 
of water, three or four digits (P,) 
Barometric pressure, inches of mercury, four 
digits (PB) 
Pressure drop across the tube bank inches of 
Mercury, four Gigi tome) 
Inlet temperature, degrees Farenheit, four 
Ga Ga ean le 
Exhaust temperature, degrees Farenheit, four 
Geqies mole) 
Tube surface temperature, degrees Farenheit, 
FOUL Gigi ts. ll 
Orifice diameter, inches (IDO) 
Speci towniinidity, Gsalnomper. pound. (ou) 
Geometry indicator ‘LGF) 
Angle indicator (LA\GLE) 


Arrangement indicator (LAF) 





1616 
1617 
1618 
1619 








1040 
150 
1670 
10670 
2090 
2200 
2250 
250 
2ef)D 
20830 
230 
Be DU 
309) 
4el() 
420 
4 e3u 
4eIl 
4e9U 
Beart) 
6020 
605i) 
10 BQ 
8 «30 
830 
Be SL 
9260 
1140 
13elu 
164006 
16010 
18240 
20050 
24200 
250 5! 
29010 
3300 
) owed) 
160 
2060 
2290 
%e5U 
4e19 
5-00 
6 «06 
6240 
ft oat 
TeZ2l 
Bel © 
830 
De AC 
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APPENDIX F 


Reduced se xpenimencal Udta 


This Appendix consists of a listing of the results of 


the data reduction. The column headings are self explanatory. 
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APPENDIX G 
CORVeCC EP OM—beocm eG diishe@mekeCrOssl On Analysis 
A regression analysis is applied to the results of the 


mata reducing computer orogram, the purpose of which is to 


fit a curve of the generalized type, Y = A + BX through 


tne reduced data. The independent variable X is the exper- 
imental value of the Reynolds Number through the test section, 
Merined throughout this text as BNRE. The independent 
Variable Y is the ratio of the Nusselt Number as measured 
experimentally to the Nusselt Number computed from the 
mcedified Grimison equation, defined by FTHETA, Tne results 
Pee chilis regression analysis will be the constants A and B, 

and the exponent M resulting in an equation relating FTHETA 
to Reynolds Number BNRE, within the range of the experimental 
Pioeas, From this equation, particularized for tne angie of 
emciination and the tube geometry, Ehe angle factcr FTHETA 
Mmey be computed, or taken from design—-oriented curves to appa 


een Sud tional tdctor te the modeled Gramiscen Eocuation, 


The analysis was Carried Out uUSing a prepared program 
available at the Applied Mathematics Laboratory of the David 
Taylor Model Basin, The program is written an Beli PCRIRAN Sa 
for the IBM 7/7090 computer. Running time is 245 minutes. The 


Program. 2S documented In Dive Repo, vue 20s7. 


Beet GeSCmp tic mae tte Bie Cr Aes et eG Chm dere (tae 


Apnpendix.e 











A maximum of 35 data points may be entered as input for 


may one Case. The program will fit a function of the form 

m= A+B es tO the data points entered by the method of least 
memares. the line of regression will be calculated and plottved@ 
Beeee Limes using three sets Gf data points, the first being 
eae entire original set of data points, the second being 

Mae set of points within two standard errors, ana the third 


Msing those points which fall within one standard error. The 


Pceandard error is defined as: 


et 


where: Y average is the average ordinate of all the points. 
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a ey tn ee ee ee 
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Yi is the generel ordinate when i is indexed from 1 to N. 


Nis the number Om spolmes. 


The method used by the program is the method of least 
Squares. An approximation tothe exponent M is made by uSing 
the first three and the last three points, averaging the 
ordinates and abcissas of these points to derive an average 
point, X average. Four points, (two greater than and two less 
than) closest to X average are taken. From these four points, 
e new average point 1s derived. The siopes of the Upper ana 


lower portion of the first approximate line of regression dre 





then determined. Dividing the lower siope by the Wpper silcper 


an expression exists from which the exponent M may be approximated. 








KR avg. = Agthy 


2 
The four values of X closest to X avg. are locaees 


x ae +X. 
iL- 


C 4-2 + X, + A + 1 
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Al 


a a + l 


in which K ee rely Cte 


Then the exponent M may be determined in 


the following manner: 
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< pmxMnt 
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Two slopes may be written from Fig. A 





mY. Y cage cake 
c 5 S 5 
eas Gl Ga BM 5 
Cc S 
\ where Yo ~ Yo is the slope at Xe + Y« 
\ am jas eae 
Xo Ke Z 
M-1 
C2 ) XY, - Ye X, + Xe 
ee a <2 
aly e 
where vy ~ Ye is the slope at Xe + X 
Xy - Ke 2 
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Bividing equation (2)—by equation (1). 


a M=-l 
(y = ve) (X-. = Xe) Xx, ae 
es os = = x 
(Xx, x) Ce Yo) oa Xe 


C 


Now let — xe Yo) (Xe = Xo) 


= D 
CX, =~ to Yo “ Yo 
Let Xr + Ko - EE 
Ko + X 
bp = Eg (ML) 


LnDb = (M-1) In E 


Then M = ine 
iigt. 0S 


which provides an approximation to M. 


+ 1 


The initial approximation of the exponent is testing for 
accuracy by considering seven valucs of the expon_nt, M being 


the center value and the others 1/4 M from this central value in 











increasing ana decreasing increments. 


HOG Gach Value ot Ltnevrexponent, the valucs of A; ana 2b 
are computed, and the sums of the squares of the derivations 
are minimized. The exponent, M, is then chosen from the set 
foes, B, and M thet has yielded the miniwum sum of the deriva-— 
Mhons squared. 


ane normal equations for theyrfunetion Y = A +4 Bx” with 


meee, ond M as the regression parameters. 


M Pa 
A ees a 
ney 8B 2X, a 
i as ee x NM LS yxy Lo 
ny a) ee 


2M 


mM 


M 
ln X, -ZY,X, Hee Sa 


SINR Ome Fs 2507, 
a i. i 
These equations are solved for A anu B by Cramer's Rule. 
The final test for the proper exponent is made by examining 
Biemcoctficlients of the three normal cqaualions, Lec. ,ecntec 
equations in two unknowns. Since A and B are unique, two of 


the three equations must be linearly dependent anu the 


meecrMinant Of ENnencocttletents wis Webec Identically eoucl tee 


ZeCLO. | 
n a Y5 
xm 2 vx," 
- a 
m x, pe a Ve ae 
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Muss PSoteue pecause the Coel ficients Of Che €qUmrion 
fmt DeLedudal wor preperitoneal) to the Goertficients Of anothems 
It can now be seen that successive iterations based on estimates 
mee unerl the determinants of Ene Coectficients are equal te 


mero to finalize the value of the coefficient M. 


The exact values of A and B can then be calculated using 


eramer's ruie. 


ine Cutput of Che program yarpeers am two forms. One fem 
is the regular output of the 7090 program which states the 
A and B coefficients and the exponent M for each of the three 
Hanes of regression. The output of a typical analysis is 
shown on pagej3,for examination. Note that with the lines of 
regression for points within two and one standard error, those 
points which were not used in the determination cf the line 
Mieee listed with the results. his list of points 1S cumulative 
mmeeugh the two cases, The adenedt ying anrormaetion which 
accompanies each gata point is placed in the computer along 
Peer the date at the choice of the User or the program. The 
Second LTOormeot output ts *through tne Sev40208Charactren. tine 
@oding for this outdbubl is “contained in arcubroucine package 
memed AMPFOS and the form of ehewoutpUuters microfilm copies 
oe Ene aeCtUdl lanes "Gfe regression. ee.) cloudie el) econ rete 
have little value from the practical standpoint, due to the 
selection by the program's Coding of rather cumbersome sGalesce 
ene plots do indicate the valiaity of the curve fit ena permis 
selection of one of the three lines as being most representative 


of the points. This output form is shown on page 145. 
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Results of Data Reduction 


Froures 27 thru 60 which fellow are plots of Nusselt 
Number versus Reynolas Number for individual geometries and 


angles. Curves for 0 = 90, 75, 60, 45 and thirty degrees are 


eiotted using experimentally derivea data, other curves for 


mre Geometry are also plotted using the Colburn and Grimisog 


equations. 
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PEE oNDL xo 


~ Widby cis Olmert eCelOlnalenes’ stance Results 
( 


Figures 51 thru 66 are plots of the computed friction 
factor versus Reynolds Number for individual geometries and 
angles of tr. The uncertainty of data for Reynolds 
Numbers lower than 2500 was such that no points were used in 


the analysis. The process of curve fittino used was the 


method of least squards. 


Peds emegudteom iin Sie 


Mhe Frictional resistance, Geta sioutrl lle to an equation 
of the form: 

BRIEF AC T=) CONST (CNR Ey 
using the method of least squares. The logarithmic equivalent 
of the equation is taken: 

Hog FRIFACs=" log CONS (4) E Oe slog NRE 
and the normal equations developed by the usual method of 
multiplying by the coefficients of the variables CONST and EXPO 


miecach case for Gach data poine sand sunning. = Sie sncrmal 


equations so developed are: 





log FRIPFAC@= (log CNRE) VExXPO een tHUMestog CONST 
(Log -FRIPAG log CNRE S-siGiog CNRE)* EXPO + (log CNRE)log COMea@ 
These normal equations can be solved simultaneously for the log 


ef -CONST cand. fxr O- 


EXPO = (THUMB) (log FRIFAC) (log CNRE)—(log CNRE) (loq FRIFAGE 
(THUMB) (log CNRE)* + (log CNRE)* 
XCON = 


(log FRIFACG)— CicaqschRe) (exeo) 
(THUMB) 


A 
Then: CONST = antilog XCON 


The parameters necessary for expression of friction factor as a 


function of Reynolds Numbers are thus developed. 
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Bereta Oondal Resistance j~malysis Computer Program Description 
Gae computer program described herein is a FORTRAN II 
program compatible with the IBM 1620 data processing system. 
mine function Of the eregeam 1S tc fit a curve whereby fricticg# 
factor is expressed as a power of Reynolds Number using the 
method of least squares. While varicble names, input, output 
and format statements are particularized to the problem at hana 
the methou is completely general and can be usea in any 
seo plication where @ leogerivulimic Least souares analysis is 
desired with only minor program modifications. [It would be 
Pessible to include thie yp procgeamein the Basic data reduction 


program to permit reduction and analysis in one computation. 


Overall COnUrOl Ozwene preOcram 1s by a DO loop andexcama, 
the number of sets of aata being analyzed. For each data Set 
enother DO loop ing@exed Soy the number er data Dolnesmin venice 
@ata set controls the -Gateularlonc we Cfnecac derarscere tne 
program first sects the summations equal to NUMP, then accents 
@s input the frict.on facror and Reynelds mumber for the fives 
Gata point. INe naturale ogari coms Cee ene siupuce Quan ll alec. 
the product of the two logarithms ana the square of the 
hogar i chin, Of Veynota. unee are computeaq. The summations 
necessary to determine the coefticients of tne normal equations 
are computed with variable names being changed as necessary 
for proper program operation, the program then loops to repeag 
the calculelLons deserved Goeve for Cac Gaca DOInNt Until aa 


data peints Pn elie tel Mave De teeelemiiteGae mine Values OF 














| ie, 


the exponent, EXPO, ana the natural logarithm of the coefficient, 
XCON Cam then be computed for the particular data set. As no 
el Logarithm f£unee on temevat table smgthie son TRAN TE labrary 
mae value of tne @oerriereme.sCelot 1s calculated by: 
CON si Poa 
when tne values of the coefficient and exponent have been 
determined an output consisting of these quantities plus an 


meentcit ying Coue, “and them oGgadrrcim Greene reoetficient as 


painted. 


The quantities which have been computed so far represent the 
perameters necessary for expressing Friction factor as a functaem 
of Reynolus Number for a particular geometry ana angle of 
pniciinations In order to correlate the values calculated wide 
the equation developed using these parameters with test results, 
values of Reynolus Number from three thousand to eight thousand 
are inserted into the equation and the corresponding friction 


factor computed. The results are printed. 


When. the celculetien@got the friction factors 15 Comoletea 
tne program loops to -pertomm, the calcuiladmloneuror Ene new. deme 


Set until the limiting gnqexers reached, 
A program Dlock diagram, 11 sting, “semple Gutpuc., and ocutpems 
summary follow. 


Computations using this program were made at C.W. Post 


College; data was prepared at webb Institute. 
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Cc FORTRAN @ikwid PROGRAM 710 CORRELATE FRICTIGNAL RESISTANCE DATA By 


8 LEAST SQUARES 
FORMAT(14) 
FORMAT ( 214} 


l 
eZ 
Be FORMAT( EO. Ga Giese | 
4 FORMAT (1LR1,2X%,5 7HRUN CODE EXPUNENT CONS TANT LOG (Gia 


© TANT) ) 
> FORMAL ClAOy Sy hr Uo g oh et Oe Os Xs 1066 } 


SFORMAT ( PHOWGxX4s pees 1 OKs 9 oo) 

10 FOURMAT(1HO,2X%,37HREYNOLDS NUMBER FRICTION FACTOR ) 
READ 1,KASES 
DO 7 MA=1,KASES 
SUMFLN=0.0 vs 
SUMREL=0.0 
SMSRG-0. 0 
SMFGRE=0.0 
READ 2,NUMB,RNCD 
DO 6 NU=1,NUMB 
READ 3,FRIFAC,CNRE 
FLG=LOGF(FRIFAC) 
RELG=LOGF(CNRE) 
FLGREL=FLG*RELG 
SRELG=RELG*¥*2 
SUMFLG=SUMFLN+FLG 
SUMFLN=SUMFLG 
SMRELG=SUMREL+RELG 
SUMREL=SMRELG 
SMSREL=SMSRG+SRELG 
SMSRG=SMSREL 
SMFLRL=SMFGRE+FLGREL 
SMFGRE=SMFLRL 

6 CONTINUE 
THUMB=NUMB 
EXPO=(THUMB*SMFGRE-SUMREL*SUMFELN) /( THUMB*SMSRG-(SUMREL*¥* 2) ) 
XCON=((SUMFLN-(SUMREL*®EXPO))/THUMB) 
CONST=2.71828%*XCON 
PRINT 4 
PRINT 5,RNCD,EXPO,CONST»XCON 
PRINT 10 
DO 8 NRE=224 
DNRE=NRE 
ENRE=DNRE*(1.0E3) 
FRICT=CONST*( ENRE**EXPO) 
PRINT 9,ENRE,FRICT 

8 CONTINUE 

7 CONTINUE 
STOP 
END 
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